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Abstract 

The rare decays B° — ► D~ K + and B + — > DfK° can occur only via annihilation 
type diagrams in the standard model. We calculate these decays in perturbative QCD 
approach. We found that the calculated branching ratio of B° — > D~ K + agreed with 
the data which had been observed in the KEK and SLAC B factories. While the decay 
B + -> D+K has a very small branching ratio at O(10 8 ), due to the suppression 
from CKM matrix elements |K* b V^|. 
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1 Introduction 



The rare B decays are being measured at B factories in KEK and SLAC. The generalized 
factorization approach has been applied to the theoretical treatment of non-leptonic B decays 
|L|. It is a great success in explaining many decay branching ratios [Q, [|. The factorization 
approach (F.A.) is a rather simple method. Some efforts have been made to improve their 
theoretical application and to understand the reason why the F.A. has gone well [||. [|. 
One of these method is the perturbative QCD approach (PQCD), where we can calculate 
the annihilation diagrams as well as the factorizable and nonfactorizable diagrams. 

The rare decays B — > D S K are pure annihilation type decays. In the usual F.A., this 
decay picture is described as B meson annihilating into vacuum and the D s and K mesons 
produced from vacuum then afterwards. To calculate this decay in the F.A., one needs 
the D s — > K form factor at very large time like momentum transfer 0{Mb\ However the 
form factor at such a large momentum transfer is not known in F.A. This makes the F.A. 
calculation of these decays unreliable. The annihilation amplitude is a phenomenological 
parameter in QCD factorization approach (QCDF)[|J, and the QCDF calculation of these 
decays is also unreliable. In this paper, we will try to use the PQCD approach, where the 
annihilation amplitude is calculable, to evaluate the B — > D S K decays. By comparing the 
predictions with the experimental data, we can test the PQCD evaluation of the annihilation 
amplitude. 

A W boson exchange causes bd — > cu or bu — > dc, and the ss quarks included in D S K 
are produced from a gluon. This gluon attaches to any one of the quarks participating in 
the W boson exchange. This is shown in Figure [I]. In the rest frame of B meson, both s 
and s quarks included in D S K have 0(Mb/2) momenta, and the gluon producing them also 
has q 2 ~ 0(M B /A). This is a hard gluon. One can perturbatively treat the process where 
the four quark operator exchanges a hard gluon with ss quark pair. It is just the picture of 
PQCD approach. 

In the next section, we explain the framework of PQCD briefly. In section we give the 
analytic formulas for the decay amplitude of B — > D S K decays. In section |], we show the 
predicted branching ratio from the analytic formulas and discuss the theoretical errors. It 
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is found that the prediction are is in good agreement with the data, and PQCD approach 
correctly gives the annihilation amplitude in consequence. Finally, we conclude this study 
in section |5|. 

2 Framework 

PQCD approach has been developed and applied in the non-leptonic B meson decays 

H for some time. In this approach, the decay amplitude is separated into soft($), hard(i7), 

and harder(C) dynamics characterized by different scales. It is conceptually written as the 

convolution, 

Amplitude- Jd 4 hd 4 k 2 d 4 k 3 Tr[C{t)$ B {k 1 )$ Da {k 2 )$K{k3)H{k 1 ,k2,k 3 ,t)], (1) 

where k^s are momenta of light quarks included in each mesons, and Tr denotes the trace 
over Dirac and color indices. C(t) is Wilson coefficient which results from the radiative 
corrections at short distance. In the above convolution, C(t) includes the harder dynamics 
at larger scale than Mb scale and describes the evolution of local 4-Fermi operators from 
mw, W boson mass, down to t ~ 0{\/KMb) scale, where A = Mb — rn^. H describes the 
four quark operator and the spectator quark connected by a hard gluon whose q 2 is on the 
order of AM#, and includes the 0(\/AMb) hard dynamics. Therefore, this hard part H can 
be perturbatively calculated. $m is the wave function which describes hadronization of the 
quark and anti-quark to the meson M . While the H depends on the processes considered, 
$az is independent of the specific processes. Determining <3>m in some other decays, we can 
make quantitative predictions here. 

We consider the B meson at rest for simplicity. It is convenient to use light-cone coordi- 
nate (p + ,p~ , p T ) to describe the meson's momenta, where p ± = -^(p ±p 3 ) and = (p 1 ,^ 2 ). 
Using these coordinates we can take the B, D s , and K mesons momenta as Pi = ^f-(l, 1, Ot), 
P 2 = ^|(i ; r 2 ,0 T ), and P 3 = %(0,1 - r 2 ,0 T ), respectively, where r = M D JM B and we 
neglect the K meson's mass M&. Putting the light (anti-)quark momenta in B, D s and K 
mesons as k\, k 2 , and k 3 , respectively, we can choose k\ = (xiP^, 0, kiy), k 2 = (x2P 2 + , 0, k 2 r), 
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and k 3 = (0,xsP 3 , k 3 r). Then, integration over k x , k 2 , and k 3 in eq.([L]) leads to 
Amplitude ~ J ^1^2^X361^6162^2^3^3 

Tr [C(t)$ B (x u 6i)$ Ds (x 2 , &2)Mz 3 , &3)^(^, 6<, t)^) e~ s(t) ] , (2) 

where 6, is the conjugate space coordinate of k iT , and t is the largest energy scale in H, 
as the function in terms of Xj and 6j. The large logarithms (In m^/i) coming from QCD 
radiative corrections to four quark operators are included in the Wilson coefficients C (t) . The 
large double logarithms (ln 2 Xj) on the longitudinal direction are summed by the threshold 
resummation||, and they lead to S t (xi) which smears the end-point singularities on Xj. The 
last term, e~ s<yt \ contains two kinds of logarithms. One of the large logarithms is due to 
the renormalization of ultra-violet divergence lnt6, the other is double logarithm In 2 b from 
the overlap of collinear and soft gluon corrections. This Sudakov form factor suppresses the 
soft dynamics effectively JTTi | . Thus it makes perturbative calculation of the hard part H 



applicable at intermediate scale, i.e., Mb scale. We calculate the H for B — > D S K decays in 
the first order in a s expansion and give the convoluted amplitudes in next section. 

3 Analytic formula 

3.1 The wave functions 

In order to calculate analytic formulas of the decay amplitude, we use the wave functions 
®M,ap decomposed in terms of spin structure. In general, $m,o/3 having Dirac indices a, j3 are 
decomposed into 16 independent components, l a/ 3, 7^, a^, (7^75)^/3, 750/3. If the considered 
meson M is B or D s meson, to be pseudo-scalar and heavy meson, the structure (7 /i 7s)a/3 
and 75 a/ 3 components remain as leading contributions. Then, $M,a/3 is written by 

$ Af,a/3 = -y== { ( pMlb)al3(l>M + TW^m} , (3) 



where iV c = 3 is color's degree of freedom, Pm is the corresponding meson's momentum, and 



b^f are Lorentz scalar wave functions. As heavy quark effective theory leads to 0^ ~ Me a 1 



then B meson's wave function can be expressed by 



®B,«p(x, b) = [(fiTsW + M B75a/ 3] <Mx, 6). (4) 

v 



According to ref . |TTJ , a pseudo-scalar meson moving fast is parameterized by Lorentz scalar 
wave functions, (ft, <p p , and 4> a as 

(D; (P)\s(z) W c(0)\0) ~ -if Ds P, [dx e ixP *<j>(x), (5) 

Jo 

(D;(P)\s(z) % c(0)\0) = -if D m 0Da [ dx e ixP *<f> p (x), (6) 

Jo 

i ( M 2 \ r 1 

(Dj(P)|^) 7 5^c(0)|0) = -f Ds m 0Ds 1 - (P„z„ - P v z») \ dx e lxPz Mx), (7) 

V m 0D s / JO 

where m £> s = ^r> s /( m c + m s ). We ignore the difference between c quark's mass and D s 
meson's mass in the perturbative calculation. This means, putting A' = Mjj a — rn c , the 
terms proportional to A'/Md s are neglected. In this approximation, the contributions of 
eq.(^|) are of higher power than those of eqs.(|, |fj) by O (j^-^ because of the factor 1 — 
in eq.(^), and we neglect the 75(7^ component in the D s meson's wave function. In addition, 
the eq.(H), eq.(|), and the relations 

jL{D;(P)\s(z) W c(0)\0) = im s (D;(P)\s(z) l5 c(0)\0), (8) 

<} (D;(P)\s(0)Wc(z)\0) = im c (D;(P)\s(0)Kc(z)\0), (9) 



dZf, 



with equations of motion lead to 

4> p {x) = 4>{x) + o . (10) 

Therefore the D s meson's wave function can be expressed by one Lorentz scalar wave func- 
tion, 

^D s M^ b ) = ~y= [(75 A)a0 + M Dslba p\ <f) Ds (x,b), (11) 



where 4>d s is defined by 



Mx) = 2jm* [x) = dk Ux) - (12) 



The wave function </>m for M = B, D s meson is normalized by its decay constant f J 



M 



f. 



d.r <> u (.r.h = 0) = (13) 



In contrast to the B and D mesons, for the K meson, being light meson, the 75a 

M 2 

component remains because the factor corresponding to 1 in eq.(0) is 0(1). Then, 

m 0D s 



K meson's wave function is parameterized by Lorentz scalar wave functions <\yjf )T as 

r dz + 

®k^(xM = P 3 / -^e- ik ^(K+(P 3 )\u pi {z)s aj (0)\0) 



75 p3<f>K(x3, 63) + m 0K -f 5 (p^(x 3 , 63) 

+ moK^/) fx - l)0^(x 3 , 63) (14) 

- a/3 

where v = (0, 1, Ot) oc P3, n = (1, 0, Ot) oc z, and mo^- = M\j{m u + m s ). In the numerical 
analysis we will use 0^-' P,T which were calculated from QCD sum rule ||11|| . They will be 
shown in section f|. 

3.2 B° -> decay 

We first consider the neutral B° decay B° — > D~K + . The effective Hamiltonian at the scale 



lower than M\y related to this decay is given as [12 



h cS = [cmom + c 2 ( /U )o 2 ( / i)] , (15) 

Oi = (bd)v-A(uc) V - A , 2 = (bc)v-A(ud) V -A, (16) 

where Ci j2 (/i) are Wilson coefficients at renormalization scale fi, and summation in SU(3) C 
color's index a and chiral projection, ^ Q ,9a7 v (l — 75)<?L> are abbreviated to (gg')y-A- The 
lowest order diagrams contributing to B° — > D~K + are drawn in Fig.^| according to this 
effective Hamiltonian. As stated above, B — > D S K decays only have annihilation diagrams. 

We get the following analytic formulas by calculating the hard part H at first order in 
a s . Together with the meson wave functions, the amplitude for the factorizable annihilation 
diagram in Fig.||(a) and (b) results in Fa~ 2 \ 

F® = -IQitC f M 2 b / dx 2 dx 3 / b 2 db 2 b 3 db 3 Ds (x 2 ,6 2 ) 
Jo Jo 

x {(1 - r 2 ) (1 - 2r 2 - (1 - r 2 )x 3 ) (f>£(x 3 , b 3 ) 
+ r (3 - r 2 - 2(1 - r 2 )x 3 ) r K (f>^(x 3 , b 3 ) 
- r(l - r 2 )(l - 2x 3 )r K </%{x 3 , b 3 )}E i f {t 1 a )h a {x 2 , x 3 , b 2 , b 3 ) 
-{(l-r 2 )x 2 ^i(x 3 ,b 3 ) 
+ 2r(l - r 2 + x 2 )r K ^(x 3 , b 3 )} E){t 2 a )h a {\ - x 3 , 1 - x 2 , b 3 , b 2 )} , (17) 
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where Cf = 4/3 is the group factor of SU(3) C gauge group, and rx = ttiok/Mb, and the 
functions E^, t\' 2 , h a are given in the appendix. The explicit form for the wave functions, 
ipMi is given in the next section. The amplitude for the nonfactorizable annihilation diagram 
in Fig.||(c) and (d) results in 

J /*1 POO 

M a = —==Q4irC F Ml / dxidx 2 dx 3 / b 1 db 1 b 2 db 2 4>b{x\, &i)0d s (x 2 , b 2 ) 
y^J^c Jo Jo 

x {(1-r 2 ) ((l-r 2 )(l-a;3)+r 2 x 2 )^(x3,62) 
+ r (x 2 + (1 - r 2 )(l - x 3 )) r K cj)^(x 3 , b 2 ) 
+ r (x 2 - (1 - r 2 )(l - x 3 )) r E -0^(x3,6 2 )}-Em(tm)/ii 1) {xi, x 2 , x 3 ,b h b 2 ) 
-{(1-r 2 ) ((l+r 2 )x 2 -r 2 ) <p A K {x 3 ,b 2 ) 
+ r(2 + x 2 + (1 - r 2 )(l - x 3 ))r K (j)^(x 3 , b 2 ) 
+ r(-x 2 + (1 -r 2 )(l ~ x 3 ))r K (t) T K (x 3 , b 2 )}E m (t 2 n )h^(xi,x 2 ,x 3 ,b 1 ,b 2 

where x\ dependence in the numerators of the hard part are neglected by the assumption 
X\ <C x 2 ,x 3 f\. The total decay amplitude for B° — > D~K + decay is given as 

A = f B F^ + M a , (19) 

where the overall factor is included in the decay width with the kinematics factor. The decay 
width is expressed as 

r(fi° - d;k+) = ^(i - r 2 )\v: b v ud A\\ (20) 

The decay width for CP conjugated mode, J3° — > DfK~, is the same value as B° — > D~K + , 
just replacing V* b V u d with V c bV* d . Since there is only one kind of CKM phase involved in the 
decay, there is no CP violation in the standard model. 



1 We don't apply this approximation to the denominators of the propagator which are sensitive to the 
variable x\. Because such a x\ behaves as the cut off, the resultant branching ratio is smaller than it given 
in ref. [T^| where the whole x\ in the hard part are neglected. 
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3.3 B+ -> D+K° decay 

The effective Hamiltonian related to B + — > D+TT decay is given as 

ffeff = ^fOk [CMOM + cmom] , (21) 

Ox = (bd) v ^ A (cu) V -A, 2 = {bu) v ^ A (cd) V -A- (22) 

The amplitude for the factorizable annihilation diagram results in — Fa~ l \ The amplitude 
for the nonfactorizable annihilation diagram results in 



i r 1 r c 

M' a = —==QA-kC f M 2 b I dxidx 2 dx 3 / b 1 db 1 b 2 db 2 4>B(xi,b 1 )(j) Ds (x 2 ,b 2 ) 
V*N C Jo Jo 



x {(1 -r 4 )x 2 (f)£(x 3 ,b 2 ) +r (x 2 + (1 -r 2 )(l - x 3 )) r K (f)^(x 3 ,b 2 ) 
+ r (-x 2 + (1 - r 2 ) (1 - a*)) r K f K {x 3 , b 2 )} E'J^Jh^ (x u x 2 , x 3 , h, b 2 ) 
- { (1 - r 2 ) ((1 - r 2 )(l - x 3 ) - r 2 + r 2 x 2 ) <^(x 3 , b 2 ) 
+ r (2 + x 2 + (1 - r 2 )(l - x 3 )) r^0£(x 3 , 6 2 ) 
+ r ( X2 _ (i _ r 2)(i _ X3 )) r rf(x 3 , 6 2 )}^(0^ 2) (^i^2, x 3 , &i, 62)] • (23) 



Thus, the total decay amplitude A f and decay width T for B + — > D+K decay is given as 

A' = -f B FU + M' a , (24) 
r(B+ - = ^(1 - r 2 )|a/| 2 (25) 



The decay width for CP conjugated mode, B — > £> s K , is the same value as _B + — > -D+if . 
Similar to the 5° decay, there is also no CP violation in this decay within standard model. 



4 Numerical evaluation 

In this section we show numerical results obtained from the previous formulas. At the 
beginning, we give the branching ratios predicted from the same parameters and wave func- 
tions that are adopted in the other works. Secondly, we show the theoretical errors due to 
uncertainty of some parameters. 
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For the B meson's wave function, there is a sharp peak at the small x region, we use 



!>s(x, b) = Nbx 2 (1 — x) 2 exp 



Mix 2 1, 



2u; b J 2 



(26) 



which is adopted in ref. [§, |J. This choice of 5 meson's wave function is almost a best fit 
from the B — > i^7r, 7r7r decays. For the Z) s meson's wave function, we assume the same form 
as .D meson's one except for the normalization [14" 



<P Ds (x) = -^=f Ds x(l ~ + " 2a:)}. (27) 

Since c quark is much heavier than s quark, this function is peaked at c quark side, i.e. 

APT 

small x region, too. The wave functions <j) K ' ' of the K meson are expanded by Gegenbauer 
polynomials, 



= ^^|l + a^-i(-l + 3a + < 2 ^(3-30^ 2 + 35a 
f K {x) = — ^=(l-2x){l + a^-3(-3 + 5e 2 )}, (30) 



where £ = 2a; — 1. The coefficients ai2pi P 2,T calculated from QCD sum rule are given in 
ref. |jl HI , and their values are 

af = 0.17, of = 0.2, = 0.212, ag = -0.148, = 0.0527, (31) 

for max = 1.6 GeV. In addition, we use the following input parameters: 

M B = 5.28 GeV, M Da = 1.969 GeV, (32) 
f B = 190 MeV, f K = 160 MeV, f Da = 241 MeV, (33) 
m 0K = 1.6 GeV, uj b = 0.4 GeV, a Ds = 0.3. (34) 

With these values and eq. (|T3"D we get the normalization factor Nb = 91.745 GeV. We show 
the decay amplitudes calculated with the above parameters at Table [3]. For the neutral decay 
B° — ► D~K + , the dominant contribution is the nonfactorizable annihilation diagrams, where 
the contribution M a is proportional to the Wilson coefficient C^t), which is of order one. The 
factorizable annihilation diagram contribution is proportional to a 2 = C\ + C2/3, which is 



one order magnitude smaller. For the charged decay B + — > DfK", it is the inverse situation. 
The Wilson coefficient in M' a is Ci(t), which is smaller than the one in Fa , cl\ = Ci/3 + C2. 
Both amplitudes A and A' are at the same order magnitude here. 

The propagators of inner quark and gluon in Figure 2 and 3 are usually proportional 
to One may suspect that these amplitudes are enhanced by the endpoint singularity 



around X{ ~ 0. This can be explicitly found in eq. (|58|, |59), where the Bessel function Yo 
diverges at X{ ~ or 1. However this is not the case in our calculation. First we introduce the 
transverse momentum of quark, such that the propagators become 1 / (xiXj + fcf.) . Secondly, 
the Sudakov form factor Exp [-S] suppresses the region of small k\. Therefore there is no 
singularity in our calculation. The dominant contribution is not from the endpoint of the 
wave function. As a proof, in our numerical calculations, for example, an expectation value of 
a s in the integration for M a results in (a s /7r) = 0.10, which gives the dominant contribution 
to B° — ► D~K + decay. Therefore, the perturbative calculations are self-consistent. 



B° 


- D-K+ 


B + 


- DfK° 


f p(2) 


-0.84 + 1.57 a 


-jB^a 


-17.16 + 6.59z 


M a 


1.99- 18.36i 


K 


9.37 + 8.34i 


A 


1.15 - 16.79z 


A' 


-7.80 + 14.93 i 



Table 1: Amplitudes(10 GeV) with parameters eqs. fl32H3l| ) 



Now we can calculate the branching ratio according to eqs. ([T^, |20|, |2^). Here we use 
CKM matrix elements and the life times ]T5|, 

(35) 
(36) 
(37) 



\V ud \ = 0.9734 ± 0.0008, \V ub \ = (3.6 ± 0.7) x 10" : 
\V cb \ = (41.2 ± 2.0) x 10" 3 , \V cd \ = 0.224 ± 0.016 
r B ± = 1.67 x 10~ 12 s, t b o = 1.54 x 10~ 12 s. 



The predicted branching ratios are 

Br(B° 
Br(B + 



> D+K°) 



4.57 x 10~ 5 , 
2.01 x 10~ 8 . 



(38) 
(39) 
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The B° — > D s K + decay is observed at Belle [16] and BaBar[IT], 



Br(5 -> D;K + ) = (4.61J;? ± 1.3) x 10~ 5 , 
Br(B° -> £>7-fr + ) = (3.2 ± 1.0 ± 1.0) x 10~ 5 , 



Belle, 
BaBar. 



(40) 
(41) 



For the B + — > DfK decay, there is only upper limit given at 90% confidence level [IS 



Br(5+ -> DjK ) < 1.1 x 10- d . 



(42) 



It is easy to see that our results are consistent with the data. 

Despite the calculated perturbative annihilation contributions, there is also hadronic 
picture for the B° -> D'K + decay: B° -> J D-7r+(p+) then J D-7r+(p+) -> D~K + through 
final state interaction. Our numerical results show that the PQCD contribution to this decay 
is already enough to account for the experimental measurement. It implies that the soft final 
state interaction is not important in the B° — ► D~ S K + decay. This is consistent with the 



argument in ref. [18[ . 

The branching ratios obtained from the analytic formulas may be sensitive to various 
parameters, such as parameters in eqs.([3~T|), eqs.(34). It is important to give the limits of the 



branching ratio when we choose the parameters to some extent. Table § shows the sensitivity 
of the branching ratio to 30% change of parameters in eqs.([H|), eqs.(Q). It is found that 
uncertainty of the predictions on PQCD is mainly due to rriQK and Ub- Below we show the 
limits of the branching ratio within the suitable ranges on vtiqk and cj&. The appropriate 



extent of itlqk can be found in ref. ITT 



1.4 GeV < m 0K < 1.8 GeV. 



(43) 



a 



K 

pl,p2,t 



in the wave functions are given as functions with respect to m^x, i an d some 



K 



input parameters, 773, uj 3 in ref. [TT[. Within eq. (pE3"|), the branching ratios normalized by the 
decay constants and the CKM matrix elements result in 

(4-57^) 



Br(B° - 
Br(B+ 



D-K 
-> D+K 



-0.26^ x 1Q - 



(2.01 



-0.16N 
-0.18; 



x 10 



-8 



190 MeV • 241 MeV J VO-0412 ■ 0.9734 / 

2 / I T T T I \ 2 



VI v ud \ 



(44) 



190 MeV • 241 MeV 



K b v cd \ 



0.0036 ■ 0.224 



(45) 
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af 


(io- 5 ) 


Br(5 H 


:io- 8 ) 

n 

- -> DfK) 


0.119 


5.01 




2.32 


0.170 


4.57 




2.01 


0.221 


4.28 




1.73 


af 


Br(5° -> £> 6 7^ + ) 


Br(5 H 


n 

- -> DfK) 


0.14 


4.38 




1.74 


0.20 


4.57 




2.01 


0.26 


4.93 




2.42 


m 0K 


Br(5 -> D~K+) 


Br(E H 


n 


1.12 


4.36 




1.61 


1.60 


4.57 




2.01 


2.08 


5.18 




2.53 






Br(5 H 


n 

- - DfK°) 


0.28 


6.36 




2.39 


0.40 


4.57 




2.01 


0.52 


3.34 




1.77 




Br(5° -> Z?jiT+) 


Br(5 H 




0.21 


4.38 




1.87 


0.30 


4.57 




2.01 


0.39 


4.95 




2.19 



Table 2: The sensitivity of the branching ratio to 30% change of parameters in eqs. 



eqs. (134]). Here we don't present the sensitivity to a5„ 2 T because the branching ratios are 
insensitive to them. 
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From the B — ► K transition form factor f+(0), the appropriate extent of u>b can be obtained. 
f+(0) calculated from PQCD at mox = 1.6 GeV in the region 

0.35 GeV < oo b < 0.46 GeV, (46) 



is consistent with f+(0) by QCD sum rules given in ref. [11]. The branching ratios calculated 
at the region of eqs. d46|) are found within 

2 / It T T I \ 2 



(47) 



Br(B ° D; ,n = (4.5735) x io- ( 190M ^.^ llfeV ) ( o^.al^ y 

Br(5+ -> = (2.01+°'?°) x 10- 8 ( 1 f |V " fe 1 . (48) 

v s J K -o.io ) ^90 MeV-241 MeVy V 0-0036 • 0.224/ v ; 

In the factorizable contribution, the B meson's wave function is integrated and normalized 
by the decay constant Thus, in the factorizable dominant decay, B + — > Dfl^ , its 
branching ratio is insensitive to the change of uJb. 



5 Conclusion 

In two-body hadronic B meson decays, the final state mesons are moving very fast, since 
each of them carry more than 2 GeV energy. There is not enough time for them to exchange 
soft gluons. The soft final state interaction is not important in the two-body B decays. This 
is consistent with the argument based on color-transparency JOJ. We thus neglect the soft 
final state interaction in the PQCD approach. The PQCD with Sudakov form factor is a self- 
consistent approach to describe the two-body B meson decays. Although the annihilation 
diagrams are suppressed comparing to other spectator diagrams, but their contributions are 
not negligible in PQCD approach || 0. 

In this paper, we calculate the B° — > D~K + and B + — > Dfl^ decays, which occur purely 
via annihilation type diagrams. The branching ratios are still sizable. The B° — > D~K + 
decay has been observed in the B factories fTJJ, 17]. This is the first channel measured in 



B decays via annihilation type diagram. The fact that the predicted branching ratio is in 
good agreement with the data means that PQCD approach gives the annihilation amplitude 
correctly, the soft final state interaction is probably small in B decays, and it is one of the 
evidences to justify PQCD approach. 
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A Some functions 

The definitions of some functions used in the text are presented in this appendix. In the 
numerical analysis we use one loop expression for strong coupling constant, 

Air 

= 1T\ 1 2/A2V ( 49 ) 

where (3$ = (33 — 2nf)/3 and nf is number of active flavor at appropriate scale. A is QCD 
scale, which we use as 250 MeV at n/ = 4. We also use leading logarithms expressions for 



Wilson coefficients Ci j2 presented in ref. [0. Then, we put m t = 170 GeV, mw = 80.2 GeV, 
and m b = 4.8 GeV. 

The function E\, E m , and E' m including Wilson coefficients are defined as 

E}(t) = a t (t)a s (t)e- S ^- S ^ t \ (50) 
E m (t) = C 2 {t)a s (t) e -SB(t)-s D (t)-s K (t) } (51) 
E' m (t) = Ci(t)a a (t) e ~SB(t)-So(t)-s K (t)^ (52) 



where 



ax(t) = ^- + C 2 (t), a 2 (t) = C x {t) + (53) 



and Sb, So, and Sk result from summing both double logarithms caused by soft gluon 
corrections and single ones due to the renormalization of ultra-violet divergence. The above 
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Sb d k are defined as 



s B (t) = 8(x 1 P+,b 1 )+2 f ^-iM, 

Ji/h ^ 

S D {t) = s{x 2 P 2 \b 3 ) + 2 f ^ 7? (/A 
S K (t) = s(x 3 P+, 6 3 ) + s((l - x 3 )P 3 + , 6 3 ) + 2 f ^7,0*0. 



(54) 
(55) 
(56) 



where s(Q,b) } so-called Sudakov factor, is given as [^0 
s(Q,b) 



Q \ ot a (fjf) 



/ Q ^f{|(27 B -l-log2)+C F log^) 

f 67 vr 2 10 2 lE \ (OLsUl'W ... 



27 







(57) 



7# = 0.57722 • • • is Euler constant, and 7 g = a s /7r is the quark anomalous dimension. 

The functions h a , and in the decay amplitudes consist of two parts: one is the 
jet function S t (xi) derived by the threshold resummationQ, the other is the propagator of 
virtual quark and gluon. They are defined by 



h a (x 2 ,x 3 , 6 2 ,6 3 ) = S t (l-x 3 ) i^-j ^ 1) (M BV /(l-r 2 )x 2 (l-x 3 )6 2 ) 
x {H£\M B ^{1 - r 2 )(l - x 3 ) b 2 )J (M B y/(l - r 2 )(l - x 3 ) b 3 )9(b 2 - b 3 ) + (6 2 <- 63)} , 



(58) 



h^\xi,x 2 ,x 3 ,b 1 ,b 2 ) 



jR { o\M B v/(1 - r 2 )a; 2 (l - x 3 ) 6 1 )J (M B ^(1 - r 2 )x 2 (l - x 3 ) b^h - b 2 ) 



+ (6l <~> b 2 ) \ X 



K (M BJ P(,-)6i; 



fOT *8)>°' 

=H«(M fl «/jf*J h), for F 2 <0 ; 



(59) 



where Hq 1 ^) = Jq (z) + iY (z), and F^s are defined by 



F 2 } = (1 - r 2 )(a; 1 - x 2 )(l - x 3 ), if 2) = x x + x 2 + (1 - r 2 )(l - x x - x 2 )(l - x 3 ). (60) 



We adopt the parametrization for S t (x) of the factorizable contributions, 

2 1+2 T(3/2 + c) 



Stix) 



VtFT(1 + c) 



[x(l-x)] c , c = 0.3, 



(61) 
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which is proposed in ref. [21|. In the nonfactorizable annihilation contributions, St(x) gives 



a very small numerical effect to the amplitude ||22||. Therefore, we drop St(x) in h a and ha ■ 



The hard scale t's in the amplitudes are taken as the largest energy scale in the H to kill 
the large logarithmic radiative corrections: 



t\ = max(M B v / (l-r 2 )(l-x 3 ), l/b 2 , 1/63), (62) 
t 2 a = max(M BV / (l-r 2 )a; 2 , l/b 2 , l/6 3 ), (63) 



4 = max(M BA /|F ( 2 ) |, M BV /(1 - r 2 )x 2 (l - x 3 ), l/&i, l/& 2 )- (64) 
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Figure Captions 



Figure 1: D S K are produced only by annihilation between b and d quark in the B meson 
from point of view of quark model. While W boson exchange causes bd — > cu, the ss quarks 
included in D S K are produced from a hard gluon. 

Figure 2: Diagrams for B° — > DjK + decay. The factorizable diagrams (a),(b) contribute to 
Fa 2 \ and nonfactorizable (c), (d) do to M a . 

Figure 3: Diagrams for B + — > D+K° decay. The factorizable diagrams (a),(b) contribute to 
Fa X \ and nonfactorizable (c), (d) do to M' a . 
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